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Abstract. We construct an individual-based kinematic model of rolling migratory 
locust swarms. The model incorporates social interactions, gravity, wind, and 
the effect of the impenetrable boundary formed by the ground. We study the 
model using numerical simulations and tools from statistical mechanics, namely 
the notion of H-stability. For a free-space swarm (no wind and gravity), as the 
number of locusts increases, it approaches a crystalline lattice of fixed density 
if it is H-stable, and in contrast becomes ever more dense if it is catastrophic. 
Numerical simulations suggest that whether or not a swarm rolls depends on the 
statistical mechanical properties of the corresponding free-space swarm. For a 
swarm that is H-stable in free space, gravity causes the group to land and form 
a crystalline lattice. Wind, in turn, smears the swarm out along the ground until 
all individuals are stationary. In contrast, for a swarm that is catastrophic in 
free space, gravity causes the group to land and form a bubble-like shape. In the 
presence of wind, the swarm migrates with a rolling motion similar to natural 
locust swarms. The rolling structure is similar to that observed by biologists, and 
includes a takeoff zone, a landing zone, and a stationary zone where grounded 
locusts can rest and feed. 

1 Introduction 

Biological swarms provide fascinating examples of natural pattern formation on short time 
scales, and on longer time scales may have significant ecological and environmental conse- 
quences [TJ |2]. The most dramatic example, arguably, is that of locusts, which cause famines 
worldwide. Of particular interest are species such as the African migratory locust Locusta mi- 
gratoria migratorioides and the desert locust Schistocerca gregaria, whose habitats together 
cover the vast majority of northern Africa, the Middle East, and southwestern Asia [3]. These 
locusts, like many others, exhibit an intriguing phase polymorphism. Individuals in the solitar- 
ious phase avoid social contact. In contrast, adult locusts in the gregarious phase form flying 
swarms. These swarms may contain up to 10 10 members, cover cross sectional areas of up to 
1000 km 2 , and travel up to 10 2 km per day for a period of days or weeks as they feed [5] causing 
devastating crop loss 0]. The mechanism for the switch to the dangerous gregarious phase is 
complex, and has been a subject of significant biological inquiry. A suite of factors recently has 
been implicated, including fractal geometry of the vegetation landscape [5] and mechanosensory 
stimulus of the locusts' back legs [Jj]. In this paper, we focus on a group of insects already in 
the gregarious phase and build a mathematical model for the destructive flying swarms. 
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Fig. 1. Schematic depiction of a rolling locust swarm, after [3|. See text for a description. The downwind 
direction is to the right. 



A migratory swarm develops from gregarious individuals in several stages. First, the insects 
may form organized bands that march along the ground; these have recently been studied in [7] . 
The marching group become a flying group through a complicated process whose details depend 
on environmental conditions including time of day, wind, sunshine, and air temperature [3JH]. In 
general, locust swarm formation begins with grounded individuals performing short, local flights 
and other movements which are uncoordinated with those of neighbors. In later stages, groups 
of locusts take shape, with each group's members sharing a common spatial orientation. Finally, 
these groups become coordinated with each other, and the swarm development culminates with 
a mass departure. As the swarm propagates, it forms a rolling motion [3J [5] which is the focus 
of this paper. 

Detailed accounts of the macroscopic swarm motion and the motion of individuals within 
the group are given in [3J |5J [HI [TO]. Here we give a brief summary. A schematic diagram is 
shown in Figure [l] The macroscopic direction of the swarm is aligned with the wind. Individual 
locusts within the swarm move in the following manner. First, flying locusts head downwind 
towards the front of the swarm. Locusts reaching the front perform a mass landing, heading 
downwards and slightly upwind until they reach the ground. These locusts remain grounded, 
resting, feeding, and possibly ovipositing, until the trailing front of the flying swarm passes 
overhead, at which point they are excited into a massive takeoff upwards and slightly upwind. 
They eventually turn to align with the wind and again fly towards the front of the swarm until 
the next landing, and so on. 

The only published mathematical models of migrating locust swarms appear in [TT]. This 
work partitions the locust population into airborne locusts and grounded locusts, and describes 
each group as a continuum density field depending on a one-dimensional spatial coordinate 
(aligned with the wind) and on time. The basic model consists of coupled, nonlinear partial 
differential equations for the two fields. These equations conserve the total number of locusts 
and account for the continual flux between the grounded and airborne populations. Several 
variations on the model are considered; these account for random motion and drift of flying 
locusts, slow motion of locusts on the ground, density dependent turning, and nonlocal social 
interactions of attractive-repulsive type (see, e.g., [12jV The questions of interest is whether 
there exist traveling band solutions, which consist of a compactly supported region of population 
density that propagates cohesively over time. The analytical and numerical studies reveal that 
no such solutions exist, and the authors conclude that their models are insufficient to describe 
migrating swarms. 

Inspired by the work of we construct a minimal model for rolling swarms with the 
goal of reproducing the macroscopic group structure and motion. The key ingredient in our 
model is the explicit inclusion of vertical structure and the barrier formed by the ground. The 
added spatial structure enables the model to support cohesive, rolling swarms akin to the ones 
seen in nature. One of our primary results is a prediction of what parameters allow for this 
rolling solution, which we accomplish by considering the statistical mechanical properties of the 
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swarm. The swarms found in our model not only roll, but have zones of landing, resting, and 
takeoff similar to the biologically observed ones. In addition to presenting a model for rolling 
swarms, our work addresses the general issue of a swarm interacting with a boundary, which to 
our knowledge has not been previously considered. 

The rest of this paper is organized as follows. In Section [2] we construct our full mathe- 
matical model for rolling locust swarms. The model accounts for social interactions, gravity, 
and advection in the direction of the wind. Section [3] examines a swarm experiencing social 
interactions in free space, in the absence of gravity and wind. This problem is equivalent to 
one well-studied in the field of statistical mechanics [T3J , and has also been examined in the 
context of biological swarming |12l 114] . We review the relevant results, which depend crucially 
on whether the parameters of the social interaction potential are chosen to be in the H-stable 
or catastrophic regime. In the H-stable regime, individuals form a crystalline group whose size 
grows as the population increases. The energy per particle of the system approaches a constant 
which we estimate numerically. In the catastrophic regime, individuals form a closely packed 
group whose size remains more or less constant. The energy per particle decreases linearly, and 
we find analytical expressions for bounds on this quantity. In Section [4] we incorporate gravity 
into the model and study it via numerical simulations. When the social interactions are in the 
H-stable regime, the locusts form an (approximately) crystalline group on the ground. On the 
other hand, in the catastrophic regime, they form a bubble- like shape with a layer of insects on 
the ground, a group in the air, and a void space in between. In Section [5] we add wind into the 
model, and see that in the catastrophic case, a rolling swarm forms, similar to those described 
in [3] . In Section [6] we discuss our model and our results vis-a-vis biological observations of 
locusts, and in Section [7] we conclude and mention some directions for future investigation. 



2 Mathematical model 

We classify our model as an individual-based, kinematic, deterministic one. By individual-based, 
we mean that we track the motion of each individual locust P^l [T¥l P^l [TBI [Tf| IT51 [TU] I2UI [^T] 
E2 E3] EU E3 E3 EH EE] rather than constructing a continuum density field [Ti l [29] 130 ] 131 ] 
EH [331 EH I3H IM1 [33 EH1 EH HOI IH1 IH1 113 HH S51 HS1 IHJ HH HS1 [3D] . By kinematic, we mean 
that the equation of motion for each locust is first order in time, so that the velocity is simply 
a function of the locust positions. This approach contrasts with the case where each locust 
follows Newton's law, so that the equation of motion is second order in time. Stated differently: 
we neglect inertial forces in our model (see the discussion in [T^] as well as the models in 
[IIl[l8l[T9jE2lE6]ES]|42]|46]|l8]). The validity of this second assumption will be discussed at 
greater length in Section [6] Finally, our model is deterministic. Though stochasticity plays an 
important role in many swarming systems E3 [5TJ [52] [53] [54] , we neglect the random motion 
of locusts. We briefly consider noise in Section [6] where we see that it does not qualitatively 
impact our results. 

We now build our model. There are N locusts in the group, and the i th locust has position 
Xj. The general model is 

x i = S i + v ff + v tt . (1) 

The (noninertial) forces acting on each locust are the social interactions Si, gravity v s , and 
downwind advection v , each of which we discuss below. The direction of swarm migration is 
strongly correlated with the direction of the wind [3] [TU] and has little macroscopic motion in 
the transverse direction. For simplicity, as in [llj . we ignore the transverse direction, and so 
our model is two-dimensional, i.e., x, = (xi, z{). Here x is a coordinate which is locally aligned 
with the main current of the wind and z is the usual vertical coordinate. The explicit inclusion 
of vertical structure in the model is a critical difference between this work and the work in |llj . 

Two locusts in isolation exert forces on each other according to basic biological principles 
of attraction and repulsion (see, e.g., EJ [TH 0H [55] ) . Repulsion operates very strongly over a 
short length scale in order to prevent collisions. Attraction is weaker, and operates over a longer 
length scale, providing the social force necessary for aggregation. For a review of evolutionary 
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aspects of attraction, including increased survival due to the selfish-herd effect, see [2]. We 
model the strength of these social forces using the function 

s(r) = Fe- r ' L -e- r . (2) 

Here, r is a distance, F describes the strength of attraction, and L is the typical attractive 
length scale. We have scaled the time and space coordinates so that the repulsive strength and 
length scale are unity. We assume that F < 1 and L > 1 so that repulsion is stronger and 
shorter-scale, and attraction in weaker and longer-scale. This is typical for social organisms 
[T2") . The social force exerted by locust j on locust i is 

S ij = s ( r ij) r ij (3) 

where = |x 3 — Xi\ is the distance between the two locusts and ry = (xj — Xj)/ry is the unit 
vector pointing from to Xj. We take the total social force on each locust in the swarm to be 
the superposition of all of the pairwise interactions, 

JY 

S i = £>i- (4) 

In the field of statistical mechanics, ^ is known as a Morse-type interaction. It has also 
been used in previous studies of swarming |12L 114] . In the biological setting, each exponential 
term arises from assuming a constant rate of transmission failure of sensory data per unit 
distance (the constant hazard function assumption). However, we do not intend for our choice 
of the Morse interaction to be taken too literally. It is a convenient phenomenological model 
which, with appropriately chosen parameters, captures the essential features of attraction and 
repulsion. 

See Figure[2]for examples of ^ with different parameters. In Figure[2ja), at short distances, 
s < and so repulsion dominates. At the vertical dotted line at r 6.2, s — and so attraction 
and repulsion balance, and the social interaction induces no movement. Following [12] . we refer 
to this distance as the comfortable distance. It is given by 

LlnF 



To the right of the dotted line at r c , s > and so attraction dominates. As r — > oo, s — > + since 
the locusts are too far away to sense each other, and hence, to interact. For this example, F = 0.5 
and L = 10. Figure [2jb) is similar, but with F = 0.25 and L = 1.5. The comfortable distance 
is r c w 4.2. Even though both examples have short range repulsion and longer-range attraction 
that both decay to zero over long distances, there is an important qualitative difference betwen 
the aggregate swarm behavior observed in these two examples. We postpone its discussion until 
the next section. 

Since our model includes vertical structure, we must consider the effect of gravity. Like all 
forces in our model, gravity operates without inertia, and hence locusts in the absence of other 
effects (social interactions and downwind advection) experience free-fall at an effective terminal 
velocity which we refer to as G. Thus, 

v 9 = -Ge z (6) 

where e z is the unit vector in the positive z direction. 

As mentioned above, locusts fly more or less in the direction of the wind. Whether the flight 
is active flight in this direction, or passive advection by the wind is irrelevant in our model (sec 
also the discussion on passive versus active flying in [8]). We include a constant drift at speed 
U, and thus the advective velocity is 

v a = Ue x (7) 
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Fig. 2. (a) and (b) Social force s in ^ exerted by one locust on another as a function of their distance 
r. The broken vertical line indicates the comfortable distance, which is r c ~ 4.2 in (a) and r c m 0.77 



in (b). (c) and (d) The corresponding potentials p in I 111, where the vertical line now indicates the 
minimum. For (a) and (c), F = 0.25 and L = 1.5, which corresponds to the H-stable regime. For (b) 
and (d), F = 0.5 and L = 10, which corresponds to the catastrophic regime. For both cases, s(r) and 
p(r) are finite as r — > 0. 



where e x is the unit vector in the positive x (downwind) direction. 
Substituting Q, (JsJ) and ^ into yields our governing equation, 



N 

£■ 

3=1 



Ge z + Ue x 



(8) 



where s is given by (2) and G, U, F and L are parameters. 

We must also define the behavior of locusts on the ground. We assume a flat, impenetrable 
ground, and thus grounded locusts may not have a vertical velocity which is negative. Fur- 
thermore, since locusts rest and feed while grounded, their motion in that state is negligible 
compared to their motion in the air. The simplest boundary condition satisfying these criteria 
is the following. For any locust on the ground, if the vertical component of the velocity as 
computed on the right side of ^ is nonpositive, then the total velocity of that locust is set 
to zero. If the vertical velocity is positive, then no adjustment is performed and the locust is 
allowed to take off. That is, 



x< set to if Zi — and i,- < 0. 



(9) 



The full model is defined by ^ and In the next section, we ignore gravity and wind 
(G = U = 0) and review what is known about "free-space" swarms. 
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3 Free-space swarm 

We ignore gravity and downwind advection, setting G = U = in Q. The simple aggregation 
model that results is 

x i = ^ S (|x i -x i |)^ X l [ (10) 

which describes locusts experiencing attractive and repulsive interactions in free space, that is, 
without other effects and without boundaries, ft will be convenient to define a potential p(r) 
which is an antiderivative of ([2]), 

p(r) = —FLeT r l L + e _r (11) 

and note that the system has a free-space energy 

1 N N 

%* = i££K^). ( 12 ) 

i=i j=i 

Equation (10 1 may then be written in gradient form, 

it = -ViEf.. (13) 



A few lines of straightforward calculation verify that (13) is equivalent to ( |10| , and that Ef s is 
an energy {dEjJdt < 0). 

Equation ( 10 1 is studied in [T5], which considers the parameter space of ([2| and asks whether 
individuals form a cohesive group and, if so, what the typical distance between organisms is. 
Certain conditions on F and L must be satisfied in order to have a cohesive group; in particular, 
L must be sufficiently large and F must be sufficiently small, so that repulsion dominates at 
short distances and attraction dominates at large distances. For the elementary case of two 
organisms, if a cohesive state exists, the distance is just the comfortable distance. For the case 
of more organisms, the so-called individual distance is shown in |12j to be shorter than the 
comfortable distance. 



The system (10) is well-known in the field of statistical mechanics [13] and a recent paper 
|14j draws connections to biological swarming. The key notion is one called H-stability which 
relates to the behavior of the group as a function of the population size N, and which depends 



|, this reduces to 
L must satisfy for 



on the shape of the interaction function s(r). For our choice of s(r) in 
conditions on the parameters F and L. We state the conditions that F anc 
H-stability at the end of this section. 

We seek equilibrium solutions by performing an numerical unconstrained nonlinear mini- 
mization of the energy Ef s in ( 12 1. Given an initial condition, the algorithm finds a single local 



minimizer of Ef s However, as one might expect, multiple minimizers actually exist. We have 
verified their existence by fixing the model parameters N, F, and L, and using a variety of 
random initial conditions. The final states achieved are not identical. However, we find that for 
the values of N we study, the different minimizers are statistically indistinguishable in that their 
total size and energy are essentially identical. The local minima seem to be tightly clustered in 
the energy landscape Ef s . 

We now discuss the two important statistical mechanical regimes, following [2]. If the 
parameters F and L are chosen in the H-stable regime, then as N increases, the spacing between 
individuals approaches a finite constant. In the parlance of |12j . the group is well-spaced. Figure 
|3^ac) shows equilibria for two swarms with N = 100 and N = 1000. The interaction parameters 
correspond to those in Figure gac). We see that organisms form a crystalline hexagonal lattice 
where the nearest-neighbor distance is approximately equal for all individuals (excluding edge 
effects). As more individuals are added to the group, the inter-organism spacing is preserved and 
the group grows to cover a larger spatial region. The energy per organism Ef s /N approaches 
a constant [13] as shown in Figure Qa). This constant is the value of Ef s /N for a perfect, 



Will be inserted by the editor 



7 





-0.6 



50 






id) 






0.6 










;. •.*,"• • v' '• : • f\ 
: •.'•* ; : •' '■ '• 


















v**"* '**•**; **• 








V. "."*.*.• .*• •. • . •"•,"*/*,** " •"•"•*.*•*"• 






-0.6 




50 









-50 







50 



-0.6 







0.6 



Fig. 3. Equilibrium free-space swarms obtained by numerically minimizing the energy Ef a given by 
(12 1. In (a) and (c), N = 100 and iV — 1000 respectively. The social interaction rule is the one depicted 
in Figure[2|ac), which corresponds to the H-stable regime. As iV is increased, the individual distance is 
preserved, so the more populous group covers a larger spatial area. In (b) and (d), the same values of 
N are used, but the interaction rule is the one in Figure J^bd), which corresponds to the catastrophic 
regime. In this iV is increased, organisms pack together more densely. 



infinite hexagonal lattice, which we have calculated numerically and indicated by the horizontal 
asymptote. 

On the other hand, if the parameters are chosen outside of the H-stable regime, the system 
is catastrophic. In this case, the energy per organism is unbounded as N — > oo, and organisms 
pack together more and more closely as the population size increases. The individual distance 
is not preserved as a function of N. Figure [3jbd) shows example equilibria, again with N = 100 
and iV = 1000; the interaction parameters correspond to Figure[2jbd). The area covered by the 
two swarms is nearly equal. As N increases, so does the density of the group. Ef s /N decreases 
linearly without bound, as shown by Figure Qb) . 

To estimate Ef s /N for the catastrophic case, we construct bounds as follows. Of the N 2 
terms inside the summation in ( 12 1, there are N self-interaction terms equal to p(0) = 1 — FL. 
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Fig. 4. (a) Energy per locust Ej s /N for equilibrium free-space swarms {G — U 

H-stable regime. The energy is given by \12\ , and the social interaction parameters F, L are those from 
Figure |}ac). As TV increases, the equilibrium configurations have a value of Ef a /N that approaches 
that of a perfect, infinite hexagonal lattice as indicated by the broken asymptote, (b) Like (a), but 
the parameters are those from Figure [2pbd) which are in the catastrophic regime. Ef s /N decreases 
linearly. Lower and upper bounds are given by ( |14[ ) and ( 15 1. These are shown (respectively) as dotted 
and broken lines, but are difficult to distinguish from the actual data on this plot, (c) Relative error in 
using the lower (dotted) and upper (broken) bounds as estimates for the data in (b). 



To find a lower bound, we may assume that the remaining TV 2 — N terms take on the smallest 
possible value, namely p(r c ), the value of the potential at the comfortable distance. Even though 
this physical configuration of locusts cannot be realized, it gives a lower bound on the energy 
per particle, namely the line 



1 



E fs /N > - p(0) + (N- l)p(r c )] = - (1 — FL) - F^(L - 1)N 



(14) 



Since the equilibria minimize the energy at least locally, and since the local minima are 
tightly clustered in energy space, almost every physically realizable configuration of locusts 
will serve as an upper bound. We construct the bound by considering a convenient state with 
N/2 locusts superposed at one point and N/2 locusts superposed at the comfortable distance 
r c . There are N self-interactions each contributing p(0). There are are additional N(N/2 — 1) 
contributions of p(0) due to interactions of locusts with others at the same site. Finally, there 
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axe N 2 /2 contributions of p{r c ). Thus, an upper bound on the energy per particle is 



Ef s /N < l - [p(0) +{\n- l)p(O) + \Np(r L )] = ±N [(1 - FL) - (L - 1) 



(15) 



The lower and upper bounds are shown (respectively) as dotted and broken lines in Figure Qb). 
Because the bounds are difficult to distinguish from the actual values of Ef s /N on the plot, 
we have plotted the relative error of the bounds in Figure [3|c) . The upper bound calculation 
is aided by the fact that the total energy of two locusts at a distance d is proportional to 
p(0) + p{d) which is negative. Our sample state is constructed by gradually superposing locusts 
on the two original ones, and this leads to the line in (15 1. If our original "base" state were 



to have positive energy, our upper bound would have been a line with positive slope, which is 
indeed an upper bound, but is not useful. In this case, one would wish to construct the bound 
differently. This will be a subject of future work. 

The biological relevance of H-stability versus catastrophe is that many organisms are thought 
to have a preferred spacing (see the review in [IS]) more or less independent of population size, 
and this corresponds to H-stability. Some other organisms are conjectured [14 to pack more 
tightly with increasing population size, which would be catastrophic behavior. The implication 
of H-stability versus catastrophe is one of the main points ultimately explored in this paper. 

The conditions on F, L in (|2| required for a cohesive and H-stable group are reviewed in 
|13j . We omit the statistical mechanical details here. The condition is 



1 < L < F 



-1/3 



(16) 



If L > F -1 / 3 > 1, then a cohesive but catastrophic group will form (if L < 1 then repulsion 
dominates at large distances and a cohesive group will not form) . Note that the condition for H- 
stability is quite restrictive. For instance, consider an attractive length scale L five times as long 
as the repulsive one. Then for the system to be H-stable, the typical force F due to attraction 
must be less than 0.8% of the typical repulsive force. We stress that one must consider the 
statistical mechanical condition (16) to correctly understand the group. The examples on the 
left and the right in Figure [2] look qualitatively similar, with short range repulsion, long range 
attraction, and a comfortable distance in between. Nonetheless, the examples in Figures [3] and 
[4] show that in the limit of large N, the systems have quite different behavior. 



4 Swarm under the influence of gravity and a lower boundary 

We expand the investigation in the previous section by considering a swarm in the presence 
of gravity (G ^ 0). If there were no boundary, the system would retain its gradient character, 
with a modified energy 

JV 

E 9 = E fs + J2 Gz i 

t=l 

and the system could be converted into the free space system by moving into a reference frame 
descending at a speed G. However, if we incorporate a horizontal boundary, the swarm descends 
and some individuals land, as shown in Figures [5] and [6] The individuals in free space still have 
velocities given by x, = —ViE y but the individuals on the ground are no longer mobile. Note 
that as an alternative, we could formally include the constraint that Zi > for grounded locusts 
by making the energy E g infinite for Zi < 0. This approach would allow locusts to move along 
the ground with the horizontal component of their free space velocity, and the dynamics would 
locally minimize the constrained energy E g . 

However, as discussed in Section|2] we choose the more realistic boundary condition ^ since 
the motion of grounded locusts is minimal. Because this "sticky" boundary condition prohibits 
horizontal motion of grounded locusts, the problem loses its gradient character. Though the 
boundary condition destroys the gradient character of the free-space problem, it is still true 
that dEg/dt < 0. The system will evolve towards minimizing E g for the locusts in the air. The 
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locusts on the ground would have a nonpositive velocity if the ground were penetrable. These 
dynamics are said to be quasi- gradient in character. 

As the system is quasi-gradient , we expect the final equilibrium state to be highly dependent 
upon the initial condition. A simple thought experiment clarifies this. Consider the case where 
all the locusts are initially randomly distributed on the ground, and consider applying an 
infinitesimal perturbation to one locust (only perturbations with positive z component are 
allowed under Eq. [9]) . The equilibrium is neutrally stable since the locust will merely fall back 
to the ground, perhaps slightly displaced in the horizontal. Stated differently, all distributions 
of locusts initially on the ground are equilibrium configurations. 

For the examples below, we consider a swarm that is transiently in free-fall with some time 
to equilibrate towards its free space configuration before individuals land. When we incorporate 
the effect of wind in the next section, the geometry of the rolling swarm state we observe appears 
to be insensitive to the choice of initial condition. 

Figure [5] shows the swarm in Figure [3ja) landing under the influence of gravity of strength 
G = 0.01. The (initially) nearly circular swarm flattens into an oblong equilibrium configuration 
on the ground. Due to the H-stable character of the problem, the swarm retains an essentially 
crystalline structure. In contrast, Figure [6] shows the swarm in Figure |3jb) landing with G = 1. 
As the group lands, a dense layer of locusts accumulate on the ground. This dense layer induces 
a force which is strongly repulsive at short ranges, which in turn leads to a void in the region 
immediately above the ground. Airborne locusts occupy a region above the void. Overall, the 
shape of the swarm is bubble-like. Figure [7] compares the final states of the two examples, from 
which the qualitative difference is apparent. 

Numerical investigations for a selected set of values of F, L and G other than those used 
in Figures [5] and [6] suggest that bubbles only form in the catastrophic case, and never in the 
H-stable case. The question of exactly what conditions enable a bubble to form is a complex 
one. It is explored in forthcoming work |56j in which we formulate a continuum description of 
the swarm and study its energy minimization analytically. 



5 Swarms under the influence of gravity, a lower boundary, and wind 

Building on the previous two sections, we now consider the full model ^ and ^ which incor- 
porates social interactions, gravity, and wind. Figure [8] shows the evolution of the H-stable state 
in Figure ma) with wind U = 0.01. As the group flies downwind to the right, locusts towards 
the front of the swarm land on the ground, and no locusts on the ground ever take off. After 
sufficient time, all locusts have landed on the ground and no further movement takes place. 

This behavior may be understood by considering the social force exerted by the swarm on 
an imaginary locust, as shown in Figure |9ja). Here we consider the initial state in Figure [7]( a) 
used for the simulation in Figure [8] The horizontal axis corresponds to the horizontal position 
x, and the vertical dashed lines indicate the (initial) support of the swarm. The vertical axis is 
Zground(x), the vertical velocity that the swarm induces on a virtual locust placed on the ground 
at horizontal location x. The dotted line represents the zero thrcshhold. In most of the support 
of the swarm, the induced vertical velocity is negative. Near the edges inside the swarm, there 
is a small region of very weak positive vertical velocity. Outside the support of the swarm, the 
vertical velocity is negligible. Thus, a grounded locust that happens to be left behind by even a 
small amount will be grounded permanently, even for very small values of G. Stated differently, 
the attractive force exerted by the swarm is too weak to induce stragglers to take off. 

Figure |9jb) is analogous, but corresponds to the initial state in Figure [7jb). In this case, the 
induced vertical velocity outside of the swarm is strong and positive. A single locust left behind 
feels a strong upward attractive force, and will take off to rejoin the swarm. Thus, we expect 
from this plot that the initial state may give rise to a cohesive, traveling swarm if gravity is 
not too strong. Simulations show that this is indeed the case. With wind U = 1, the bubble- 
like state in Figure mb) undergoes a rolling motion to the right. Locusts at the front of the 
airborne swarm land on the ground, while those at the back of the grounded group take flight. 



Two snapshots of the rolling motion are shown in Figure 10 ^ab); note the qualitative similarity 
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Fig. 5. Snapshots of the H-stable swarm in Figure |3| a) landing on the ground under the influence of 
gravity. The broken line guides the eye to the top of the swarm. As the swarm lands, the macroscopic 
behavior is to (approximately) maintain a crystalline structure. Here, N = 100, the social interaction 
rule is the one in Figure J^ac) , and the gravity parameter is G = 0.01. 



to Figure [T] To further understand this motion, we examine time series of Xi and Zj for a single 
locust, shown in Figure 10 'cd). The locust flies upwards and downwind, arcing over the void 
section in the middle of the swarm, with Xj, increasing (approximately) linearly in time. The 
locust reaches a maximum height and then begins descending, with Xi still increasing linearly. 
The locust settles on the ground in the landing zone at the front of the swarm, and at this point 
its horizontal motion stops due to the boundary condition this is seen as the horizontal 
plateaux in Figure 10 c). After a time in the stationary zone on the ground, the void and the 
swarm have passed overhead, and the locust is in the takeoff zone. Nearly left behind by the 
flying swarm, it leaves the ground in order to catch up. 

One question of interest is whether all locusts undergo motion similar to that described 
above and shown in Figure [To|cd) . To answer this question, we have computed basic statistics 
describing the landing/takeoff pattern of the swarm members. This pattern is visualized in 
Figure [TT] The plot depicts 20 locusts randomly selected from the population. The vertical axis 
indexes the locusts, and the horizontal axis is time. A darkened pixel represents a time at which 
a given locust was on the ground. For the entire hundred-locust swarm, individuals averaged 
over the course of the simulation 6.7 segments on ground with standard deviation of 1.1. The 
mean duration of each grounded segment was 3.0 time units with standard deviation of 0.23. 
These figures indicate that all locusts follow a similar landing/takeoff pattern. Furthermore, 
across all times, the average fraction of the population grounded at a given time was 0.39 with 
a standard deviation of 0.0088. Across all individuals, the average fraction of time spent on 
the ground was 0.39 (since averaging over locusts commutes with averaging over time) with a 
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Fig. 6. Snapshots of the catastrophic swarm in Figure|3|b) landing on the ground under the influence 
of gravity. The broken line guides the eye to the top of the swarm. The equilibrium state consists of a 
group of locusts in the air, a dense group of locusts on the ground, and a void separating them. Here, 
N = 100, the social interaction rule is the one in Figure [2jbd), and the gravity parameter is G = 1. 



standard deviation was 0.049. These statistics confirm that in terms of distribution between 
the ground and the air, the swarm looks the same at all times. 



6 Discussion 



Several aspects of our model and our results bear further discussion vis-a-vis biological observa- 
tion. First, as discussed in Section|2] we have used the framework of a kinematic (as opposed to 
dynamic) model, and thus we have neglected inertial forces acting on locusts. While this may be 
a good approximation during the bulk of the flying stage, it is likely not accurate during takeoff 
and landing. Relatedly, takeoff and landing in our simulation occur in the upwind direction, as 
seen in Figure [10] However, in the schematic depiction in Figure [l] locusts take off and land 
in the downwind direction. This phenomenon is described in more depth in [3J, which cites a 
combination of behavioral, physiological and aerodynamic factors as explanation. It is possible 
that incorporating these factors within a dynamic model formulation would lead to simulated 
swarms with the correct landing and takeoff direction. 

Second, the social interactions in our model are of isotropic form, i.e., a given locust senses 
equally well in every direction, and the influence of any other locust only depends on the distance 
(and not the angle) between them. Sensing within a swarm is thought to be both visual and 
auditory fSJ. Isotropy is a likely a good model for hearing. For sight, we think of it as being a 
first approximation. Anisotropic interactions have been modeled in a one-dimensional setting 



Will be inserted by the editor 13 



20- 



10- 



(a)- 



• • • 



0^ . . • • . . • • • • * • ' • • • ■ ■ -L- 

-40 -20 20 40 




-1.5 -1 -0.5 0.5 1 1.5 

Fig. 7. Blow-up showing final states from (a) Figure [5] and (b) Figure [fi] Within each sub-plot, the 
scales on the x and y axes are equal. The locusts land on the ground due to gravity, but the qualitative 
character of the final state depends on the statistical mechanical regime in which the parameters F, L in 
|2| lie. In the H-stable case (a) the group has an oblong crystalline structure, while in the catastrophic 
case (b) it has a bubble-like structure with a void above the ground. 



in |49j . The effect of anisotropic interactions on locust swarming may be considered in future 
work. 

Third, our model does not incorporate random motion of locusts. In reality, flying locusts will 
experience random motion due to turbulent air currents and due to the inherent imprecision 
of their sensing and movement. As a preliminary test of the effect of randomness, we have 
conducted simulations that include noise added to the right hand side of d8J). For the i th locust, 
we chose the noise to have direction uniformly distributed from [0, 2ir\ and the magnitude 
uniformly distributed from [0, |vj|] where Vj is the deterministic portion of the velocity, i.e., 
the right hand side of ^ . Even this large amount of noise did not qualitatively change any of 
our results about rolling (or dissipated) swarms. 

Finally, our model produces rolling swarms only when the corresponding free-space group 
has a social interaction potential in the catastrophic regime. As discussed in Section [2] there is 
a widely-held belief that every species has a comfortable distance that is independent of group 
size and movement (or lack thereof) and depends only on environmental conditions; see the 
extended review in |12j . On the other hand, estimated locust swarm densities vary over three 
to four orders of magnitude |10j . Some observations found such variation over different parts 
of a given swarm, and within a given swarm over a period of a few hours, as summarized in 
[3]. These variations in density are not well-modeled by an H-stable potential. Of course, by 
using a catastrophic potential, one obtains a situation where the density of the group becomes 
infinite as N — > 00 which is obviously unbiological. Nonetheless, we conjecture that within a 
range of N, a catastrophic potential is a reasonable modeling assumption. A similar conjecture 
has been made for swarming Myxococcus xanthus cells [14] . 
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Fig. 8. Evolution ol the H-stable state in Figure^a) with wind U = 0.01. Locusts fly downwind, and 
those a the front of the group land on the ground, never to take off again. The simulation reaches an 
equilibrium once the entire group has been smeared out along the ground. 



7 Conclusions 



In this paper, we have constructed a discrete kinematic model for rolling locust swarms incor- 
porating social interactions, gravity, wind, and the effect of an impenetrable boundary, namely 
the ground. We have studied the model using numerical simulations and tools from statistical 
mechanics, namely the notion of H-stability. Our simulations suggest that whether or not a 
swarm rolls in the presence of wind depends on the statistical mechanical properties of the 
corresponding free-space swarm. For a swarm that is H-stable in free space, gravity causes the 
group to land and form a crystalline lattice. Wind, in turn, smears the swarm out along the 
ground until all individuals are stationary. In contrast, for a swarm that is catastrophic in 
free space, gravity causes the group to land and form a bubble-like structure. In the presence 
of wind, the swarm migrates with a rolling motion similar to natural locust swarms observed 
by biologists. In the rolling swarm, all individuals land approximately the same number of 
times, and spend approximately the same amount of time on the ground during each landing. 
This captures an important feature of the real locust swarms, namely that each individual has 
adequate time to rest and feed. 

Our study focused on two sets of parameters, one from each statistical mechanical regime. 
We have also conducted simulations with a limited set of other parameters and the same 
qualitative result was seen: swarms with a bubble structure roll, and those with a crystalline 
structure do not. However, this is not a proof. As mentioned in Section [4j to further understand 
when a bubble forms, we are investigating a continuum analogue of our model without wind; 
results will be published elsewhere. The development of continuum models will also be useful 
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Fig. 9. Vertical component of velocity z groun d induced by a swarm on a"virtual locust" positioned 
on the ground at location x. Where z groun d < 0, virtual locusts cannot take flight for any values of 
the gravity parameter g. Where Zground > 0, virtual locusts can take off if g < z gr0 und- The dotted 
horizontal line guides the eye to z groun d = 0. (a) H-stable case corresponding to the swarm in Figure 
[7ja) . Virtual locusts outside the swarm boundaries can take off only for miniscule g. The broken vertical 
lines indicate the left and right boundaries of the swarm, (b) Catastrophic case corresponding to the 
swarm in Figure [7^b) . The swarm induces a large upward velocity on virtual locusts outside of the 
swarm, suggesting why locusts are not left behind when the swarm migrates. 



because it will enable the analytical and computational study of larger swarms, as opposed to 
the relatively small groups we have considered here in the framework of a discrete model. 

We have presented a biologically-motivated model that reproduces many features of locusts 
swarms observed in nature. We conclude by pointing out that the parameters in this model are 
chosen heuristically. In particular, the functional form of the social interactions, while plausible, 
is at this point speculative. Our hope is that this work will inspire biologists to work toward 
gathering data that can inform and improve this model, so that eventually more quantitative 
comparisons with nature will be possible. 
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Fig. 10. Evolution of the catastrophic state in Figure[7|b) with wind (7 = 1. The group rolls downwind 
with locusts at the front of the airborne swarm landing on the ground and those at the back of the 
grounded group taking flight, as shown in snapshots at times (a) t = 5 and (b) t = 15. Note the 
qualitative similarity to Figure [I] The x and z coordinates of a single locust are tracked in (c) and (d) 
respectively. When the locust is grounded (z — 0) there is no horizontal motion due to the boundary 
condition |9|. When the locust is flying, x increases approximately linearly. 
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